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A  method  is  proposed  for  separating  measured  values  of  total 

scattering  cross  sections  for  reactive  systems  into  their  component 

elastic  and  reactive  parts.   For  scattering  of  thermal  energy  beams  the 

analysis  yields  the  long  range  Van  der  Waals  potential  constant  C 

and  the  magnitude  of  the  reactive  cross  section  Q~  .   Because  the  method 

requires  unusual  accuracy  in  the  total  cross  section,  it  is  necessary  to 

correct  the  raw  data  for  nonunif  orniity  of  the  scattering  gas  number 

density  along  the  beam  path,   This  is  accomplished  by  a  kinetic  molecular 

theoretical  calculation  of  the  steady- -state  distribution  of  scattering 

gas  for  a  scattering  chamber  with  rectangular  channel  slits  with  finite 

length.   The  correction  is  applied  as  a  factor  d  to  the  conventional 

Beer-Lambert  expression  for  scattering  in  a  homogeneous  medium.   Trie 

scheme  has  been  tested  by  collecting  and  analysing  data  for  two  systems; 

one  nor_reactive,  Cs  -  Ar,  and  one  reactive,  Cs  -*CHLC1CH„C1.   The  direct 


results  are  C(6)(Cs  -  Ar)  -  5-^8x10  5&   ere  cm  ;  G(  J(Cs  -  G^ClCr^Cl)  = 
2.3^x10    ere-  cm  ;  and  0   =  U5A2.  This  experimental  value  of  Q 
is  used  to  compute  an  upper  limit  to  the  reaction  rate  constant 
(k  =  2x10   cm3  mole"  sec"  )  and  to  estimate  the  electron  affinity  of 
CH  C1GH2C1  (E.A.  =  1.3  ev). 


viii 


CHAPTER  1 
INTRODUCTION 

In  the  early  development  of  chemical  kinetics  the  concept  of 
reactive  cross  sections  was  developed  as  a  means  of  explaining  why  not 
all  collisions  between  chemically  reactive  molecular  species  lead  to 
reaction.   Data  from  gas  phase  kinetic  studies  used  in  simple  collision 
theories  resulted  in  typically  small  (ca.  a  few  square  angstroms) 
effective  reactive  cross  sections.   This  indicated  that  the  reactants 
would  have  to  approach  each  other  within  an  angstrom  or  so  before 
reaction  could  occur.   In  early  atomic  flame  experiments  Polanyi  showed 
that  for  some  fast,  highly  exothermic  chemical  reactions  such  as 

K  +  Br2  =  K3r  +  Br, 

the  reactive  cross  sections  were  quite  high  (on  the  order  of  a  few 

1 
hundred  square  angstroms). 

The  most  successful  attempts  to  gain  better  understanding  of  how 
molecular  behavior  influences  chemical  reaction  rates  have  been  carried 
out  in  the  science  of  molecular  reaction  dynamics.  The  types  of 
experiment  which  provide  the  most  direct  measurements  of  collision 
dynamics  are  those  involving  the  scattering  of  molecular  beams. 

In  the  past  several  years  many  researchers  have  devoted  a  consid- 
erable amount  of  attention  to  the  determination  of  reactive  cross 
sections  with  the  use  of  two  intersecting  molecular  beams.   These  reactive 
cross  sections  are  extracted  from  angular  dependent,  or  differential, 


reactive  crocs  sections  IR(Q,p,U)   which,  vfhen  obtained  in  this  fashion, 
are  defined  by 

iEafi,i>)  =■ -- — .  (i) 

nA  fA(VA)   n3  f3(^)yrdA 
where  nA,  n^  and  n  are  particle  densities  in  the  incident  crossed  beams, 

AH  C 

f .  is  the  velocity  distribution  function  for  particle  A  having  velocity 

V.,  Tis  the  beam  intersection  volume,  and  V   is  the  relative  velocity. 

0  and  0  are  the  nolar  angles  which  identify  the  position  of  a  detector 

2 
of  solid  angular  width  dA.   Integration  of  the  differential  reactive 

cross  section  over  all  angles  about  the  beam  intersection  gives  the 

total  reactive  cross  section 


•2tfrtf 


QpX(W=f27    IR(6,^)  sin9d9dp    .  (2) 


Part  of  the  interest  in  reactive  cross  sections  lies  in  the  fact 
that  the  reaction  rate  constant  can  be  derived  from  the  reactive  cross 
section, .  For  the  process 

A  +  B  =  C  +  D 

the  rate  constant  k  is  defined  as 

1     dnp       1      dn 


A 


(nA)  (n  )  dt     (nA)  (n^  at 


(3) 


To  obtain  the  rate  constant  in  terms  of  reactive  cross  section  all  that 

3 
is  needed  is  a  substitution  for  dnc  from  equations  1  and  2.   Thus, 


If  all  channels  of  interaction  except  for  chemical  reaction  are  in 
Boltzmann  equilibrium,  then  the  temperature  dependent  rate  constant  k(T) 
is  computed  by 

/  -  \3/2r        \-*a    2 

k(T)   =    V%(V)   exp WZ  dV,  (5) 

\  2Tfk  t'         I  l2  k  T> 

where  u  is  the  reduced  mass  of  the  reactants. 

Determination  of  total  reactive  cross  sections  from  differential 

cross  section  measurements  is  hindered  by  several  problems  which  result 

in  only  an  estimate  of  the  order  of  magnitude  of  the  total  reactive 

cross  section.   Equation  2  shows  that  the  differential  cross  section 

must  be  determined  through  all  angles  for  calculation  of  Q  .   Designs 

of  present  apparatus  have  made  this  impossible.  Also,  since  very  little 

of  the  beam  actually  undergoes  reaction,  there  is  very  feeble  intensity 

of  the  scattered  products.   Thus,  angular  and  velocity  resolution  are 

5 
usually  sacrificed  for  a  detectable  signal. 

It  is  evident  that  the  total  cross  section  has  within  it  information 

about  all  the  scattering  processes  that  have  occured.   These  include 

elastic,  reactive  and  other  inelastic  events.   It  has  been  shown  that 

this  can  be  expressed  by  the  following  relationship 

<W<w  -  <fe<w  +  VW  +  V"  •  (6) 

where  Q~  ,  is  the  total  cross  section,  0   the  reactive  cross  section, 

Qnn  the  elastic  cross  section,  and  Q-  the  cross  section  for  the  remaining 

6 
inelastic  events.   The  total  cross  section  is  obtained  from  experiment 

by  measuring  the  attenuation  of  beam  intensity  for  a  given  scattering 
gas  number  density  in  a  scattering  box.   Precision  of  these  measurements 


is  now  better  than  1%.      The  reason  for  such  precision  is  that  intensity 
measurements  are  made  directly  in  line  with  the  beam,  which  results  in 
ample  signal  for  very  good  angular  and  velocity  resolution.  Q^Qt   is 
obtained  using  a  modified  Beer's  law  relationship 

In  (I0  /  I)  R 
«n0(L+T)FB 

where  I  and  IQ   are  intensities  of  scattered  and  unscattered  beam,  nQ 
is  the  scattering  gas  number  density,  L  the  length  of  the  scattering 

box,  and  T  the  thickness  of  the  scattering  box  slits.   R  is  an  apparatus 

7 
resolution  correction  factor,  F..   is  the  Berkling  correction  factor  for 

8  B 
motion  of  the  scattering  gas,  and  d  is  the  correction  factor  for  non- 
uniformity  in  the  scattering  gas  distribution. 

landau  and  Lifshitz  have  shown  that  for  a  potential  of  the  form 

V(r)  =  c(s)  /  rs   , 

where  C^  is  a  constant  and  r  is  the  distance  between  the  two  colliding 
particles,  the  elastic  cross  section  can  be  approximated  by  the  expression 


(s)    ,*      2/(s-D 
^QJ)   =  8.0826  (C[S)  /hV) 


(8) 


9 
where  V  is  the  relative  velocity.   Numerous  experiments  have  shown 

10 
this  to  be  a  valid  approximation  and  that  the  value  of  s  is  usually  6. 

Therefore,  it  seems,  by  obtaining  QTot  (V)   from  experiment  and 

equation  ?  and  by  calculating  Q.mQJ)   with  equation  8,  the  combination 

of  the  reactive  and  nonreactive  inelastic  cross  sections  can  be  derived 

from  equation  6.   Since  in  this  experiment  only  thermal  energy  cesium 

beams  were  used,  there  is  not  sufficient  energy  to  produce  rotational 

or  vibrational  excitations  as  would  be  needed  for  nonreactive  inelastic 


'J 

processes.   Thus,  Q^     accounts  for  most  of  this  inelastic  term. 

This  study  presents  the  experimental  apparatus  and  the  first 
results  obtained  using  this  method.   The  system  investigated  in  these 
initial  experiments  is 


Gs  +  CH2G1CH2C1  =  GsGl  +  CHg&CHg 


CHAPTER  II 
EXPERIMENTAL 

Existing  Apparatus 

The  apparatus  used  for  this  experiment  is  described  in  detail 

11,12,13.14 

elsewhere  in  the  literature.  However,  a  brief  description 

is  presented  here  for  clarification  of  the  experimental  technique. 
Figure  1  shows  the  physical  arrangement  of  the  molecular  beam  apparatus. 
The  bulkhead  located  between  the  beam  source  chamber  (oven  and  scattering 
box)  and  the  detection  chamber  (velocity  selector  and  detector)  allows 
for  differential  pumping  between  these  two  chambers.  Both  chambers  are 
pumped  with  oil  diffusion  pumps  backed  by  mechanical  vacuum  pumps,  and 
the  beam  chamber  also  has  cryogenic  pumping  with  liquid  nitrogen  traps. 

The  cesium  beam  effuses  through  the  slit  (0.005  by  O.635  cm)  of  a 
copper  oven  which  is  heated  by  several  cartridge  type  heaters.  The 
design  is  such  that  the  oven  is  loaded  with  cesium  in  a  glass  ampule 
prior  to  assembly  of  the  beam  chamber,  and  the  cesium  is  released  by 
breaking  the  ampule  when  the  apparatus  is  under  vacuum  and  ready  for 
the  experiment. 

Collimation  of  the  beam  is  provided  by  the  scattering  box  slits 
(0.010  by  O.635  cm).   As  shown  in  Figure  2,  the  slits  are  cut  in  the 
faces  of  hollow  plugs  which  can  be  moved  in  or  out  to  change  the 
scattering  path  length.   Pressure  of  the  scattering  gas  is  monitored 
in  the  cavity  below  the  slit-plugs  by  a  thermocouple  gauge. 

The  velocity  selector  is  designed  to  select  simultaneously  two 


Figure  1.  Sectional  view  of  basic  apparatus. 
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11 

beams  at  right  angles  and  in  phase.   However,  the  present  experiment 
requires  the  use  of  but  one  beam.  As  shovrn  in  Figure  1,  the  beam  is 
aimed  at  the  ed~e  of  the  horizontal  selector.  The  transmitted  portion 
reaches  the  detector  via  passage  through  spiral  grooves  cut  in  the 
bottom  face  of  the  selector  disc. 

The  detector,  a  surface  ionization  type,  is  located  in  the  hollow 
cavity  in  the  bottom  of  the  velocity  selector  disc.   Figure  3  shows  an 
isometric  view  of  the  detector  with  the  shield  removed.   The  cesium 
beam  atoms  pass  through  the  detector  to  the  hot-wire  filament  (Platinum, 
0.01  cm  diameter)  to  which  they  transfer  an  electron  and  evaporate  as 
Cs+  ions.  A  positive  potential  (with  respect  to  the  guard)  on  the 
repeller  forces  these  ions  back  through  the  slit  in  the  guard.   A 
potential  across  the  focussing  plates  causes  the  ions  to  swing  in  an 
arc  and,  thus  strike  the  collector  plate  which  is  connected  to  the 
input  of  a  Gary  Model  31  Vibrating  Reed  Electrometer. 

Modifications 

The  ultimate  success  of  this  type  of  experiment  depends  heavily 
on  the  precision  of  beam  intensity  and  scattering  gas  pressure  measure- 
ments. The  following  apparatus  modifications  were  made  to  improve  the 
reliability  of  these  measurements.   Host  of  these  improvements  are 
aimed  at  decreasing  the  time  between  measurements  of  scattered  and 
unscattered  beam  intensities  and  the  corresponding  scattering  gas 
pressures. 

Detection  System 

The  beam  intensity  detection  scheme  employed  involves  measuring 
the  time  required  for  the  electrometer  to  be  charged  between  two  fixed 


Figure  3.      Ionization  detector  with  shield  removed.      1)   Shield,    2)  collector 
(connected  to  electrometer  input),    3)  focusing  plates,   1+)  guard, 
5)   Pt-vrire  filament,    6}  filament  supports,    7)~repeller. 
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potentials.   Originally  these  potentials  were  set  by  a  series  of  micro- 
5vjitcb.es  installed  at  various  positions  along  the  slideuire  of  a  strip- 
chart  recording  potentiometer.  The  charging  time  was  measured  by  a 
universal  counter  gated  to  start  and  stop  on  signals  from  the  micro- 
switches,  and  the  total  lapsed  time  was  printed  on  a  paper  tape.   The 
new  system  works  in  essentially  the  same  manner,  but  the  mechanical 
switching  system  has  been  replaced  by  a  solid-state  electronic  circuit. 
This  device  monitors  the  output  of  the  electrometer  and  starts  and  stops 
the  counter  at  preset,  extremely  stable  potentials  set  on  voltage 
comparators.  The  device  also  automatically  grounds  the  electrometer 
input  until  the  potential  falls  below  a  certain  level,  at  which  time  the 
electrometer  input  is  automatically  unshorted  and  the  counter  reset  in 
preparation  for  the  next  intensity  measurement. 


?eam  Oven 


As  was  mentioned  earlier,  the  beam  source  oven  is  heated  by  several 
electric  cartridge  heaters.   Most  of  these  heaters  operate  at  a  constant 
current  adjusted  to  bring  the  oven  temperature  to  a  degree  or  so  below 
the  desired  temperature.   Current  through  the  remaining  heaters  is 
controlled  by  a  pyrometer  so  as  to  briag  the  temperature  up  to  that 
desired.   The  original  pyrometer  (API  Model  450-PR)  was  equipped  with 
a  mechanical  switching  arrangement  for  regulation  of  temperature.  This 
was  replaced  by  a  new  design  pyrometer  (API  Model  901-3)  which  incor- 
porates a  -ohoto-electric  system  for  temperature  regulation.   The  new 
unit  proportions  the  current-on-time  as  the  temperature  approaches  the 
desired  temperature.   This  resulted  in  improved  beam  oven  temperature 
regulation  which  in  turn  improved  the  stability  of  the  beam  intensity. 


15 

Scattering  Gas  Handling  System 

The  experimental  procedure  used  in  this  experiment  requires 
measuring  alternately  scattered  and  unscattercd  beam  intensity.  So 
that  the  time  required  to  fill  and  evacuate  the  scattering  box  could 
be  minimized,  a  remote  operated  "pump-out"  value  was  installed 
immediately  before  the  gas  input  port  to  the  scattering  box.   Figure 
4(a)  shows  an  exploded  view  of  the  valve.  Stability  of  the  valve  in 
either  the  fill  or  evacuate  mode  is  provided  by  coupling  of  the  magnetic 
field  of  the  bar  magnet  with  the  soft-iron  core  of  the  electromagnet. 
Figure  4(b)  shows  a  cross  sectional  view  of  the  valve  in  the  fill  mode, 
while  Figure  4(c)  shows  the  valve  in  the  pump-out  mode.   Installation 
of  the  valve  reduced  the  delay  time  between  the  two  modes  of  operation 
from  about  3  minutes  to  "}Q   seconds. 

In  order  to  increase  the  stability  of  the  scattering  gas  delivery 
pressure,  and,  thus,  its  number  density  along  the  beam  path,  an  insulated 
gas  handling  system  was  fabricated;  the  schematic  of  which  is  shown  in 
Figure  5.  This  eliminated  fluctuations  in  the  scattering  gas  number 
density  due  to  changes  in  the  laboratory  room  temperature. 

Experimental  Procedure 

Prior  to  each  experimental  run  the  apparatus  is  allowed  to  reach 
a  steady  condition  with  all  pumps ,. traps ,  heaters,  and  electronic 
equipment  in  the  laboratory  operating.  This  usually  requires  a  minimum 
of  48  hours.   Upon  attainment  of  this  steady  condition  the  velocity 
selector  is  turned  on  and  synchronized  at  the  velocity  corresponding 
to  the  most  probable  velocity  of  the  beam  atoms.  At  this  velocity  the 
cesium  ampule  is  broken,  thus,  initiating  the  beam.   When  the  beam  has 


Figure  4.   (a)  Exploded  view  of  pump-out  valve.   1)  Cover,  2)  bearings, 
3)  valve  core,  4)  valve  body,  5)  electromagnet  support,  6) 
bar  magnet  and  collar,  ?)  electromagnet  with  soft-iron  core, 
(b)  Schematic  of  valve  in  fill  mode.   (c)  Schematic  of  valve 
in  pump-out  mode.   Arrows  represent  gas-flow  direction. 
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Figure  5*  Schematic  of  insulated  gas  handling  system.  1)  Precision 
leak  valve,  2)  bellows  type  vibration  eliminator,  3)  hand 
expansion  valves,  4)  scattering  gas  supply  tank. 
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stabilized,  the  scattering  gas  is  introduced  into  the  scattering  "box 
and  adjusted  for  a  beam  attenuation  of  20  to  25?j.  With  this  level  of 
attenuation  at  this  velocity  there  will  still  be  adequate  attenuation 
at  the  high  velocities  and  not  overwhelming  attenuation  at  the  low 
velocities. 

Following  these  initial  steps  and  acceleration  of  the  velocity 
selector  to  the  maximum  speed  desired  for  investigation,  the  actual 
data  collection  begins.   With  the  velocity  selector  coasting  down, 
the  intensity  (i.e.  the  time  for  charging  the  electrometer  between 
two  preset  potentials)  of  the  beam  is  determined  several  times  along 
with  the  background  pressure  (eraf  of  the  thermocouple  gauge)  in  the 
scattering  box.  Gas  is  then  directed  into  the  scattering  box  and  several 
scattered  intensity  measurements  are  made  together  with  the  pressure 
in  the  scattering  box.   During  each  intensity  measurement  the  average 

angular  velocity  of  the  velocity  selector  is  being  monitored  and 

13 
recorded  automatically.    Ibis  procedure  of  measuring  series  of  data 

for  scattered  and  unscattered  conditions  is  maintained  while  the  velocity 

selector  coasts  down  through  the  desired  range  of  velocities.   The 

pressure  data  that  are  taken  are  in  the  form  of  output  eraf  of  the 

thermocouple  gauge  as  determined  by  a  Leeds  and  Northrup  Type  K-5 

Potentiometer  with  a  DC  Null  Detector  (Leeds  and  Northrup).  As  mentioned 

before,  the  beam  intensity  is  inversely  proportional  to  the  amount  of 

time  needed  for  the  beam  signal  to  charge  the  electrometer  from  one 

preset  potential  to  another  higher  preset  potential. 

Table  1  lists  typical  parameters,  times,  and  conditions  which  are 

experienced  during  an  actual  experiment.  These  data,  together  with 

those  described  above,  are  then  processed  with  extensive  aid  of  a  high 

speed  computer. 
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Table  1 
Typical  Operating  Parameters 


fnr  &rr 


Temperatures 

Oven  135  C 

Scattering  chamber  35  G 

Laboratory  25  0 


Pressures 


-7 


Background        .  ■  5  x  10"'  torr 

-4 

Scattering  gas  5  *  10   torr 

Times 

Electrometer  charging  5  -  15  sec 

Delay  between  scattered  and  unscattered  modes  3°  sec 

Total  time  elapsed  for  experiment  2.5  -  3  hrs 

Velocity  Ranges 

Selector  °000  -  1^00  rpm 

4      -1 

Bean  atoms  1  -  6  x  10  cm  sec 

4      -1 

Scattering  gas  2  x  10  cm  sec 

Dimensions 

Oven  slit  0.005  x  0.6 35  cm 

Scattering  box  slits  0.010  x  O.635  cm 

Velocity  selector  grooves  0.159  x  O.635  cm 

Detector  wire  diameter  0.010  cm 

Distances 

Oven  slit  to  facing  scattering  box  slit  12.50  cm 

Detector  to  facing  scattering  box  slit  31.05  cm 

Between  scattering  box  plugs  0.1176  cm 


22 
Table  I  -  continued 


Number  Density  Correction  Factor  OC  -  O.6775 

Apparatus  Resolution  Angle  (Kush  Angle)  5.95  x  10   radians 

Attenuations  20  -  25  % 


chapter  iii 
calculations 

Procedure 

The  overall  process  of  reducing  the  intensity,  pressure  and 
velocity  data  to  the  absolute  total  cross  section  fallowed  by  separation 
of  the  elastic  and  reactive  parts  is  quite  lengthy.   An  outline  of  the 
computational  scheme  is  presented  first,  followed  by  more  detailed 
descriptions  of  each  step. 

First  the  electrometer  charging  times  and  the  angular  velocities 
of  the  velocity  selector  are  converted  respectively  to  beam  intensities 
and  beam  velocities.   Using  an  arbitrary  function  of  the  bean  velocity 
the  beam  intensities  are  reduced  in  order  to  flatten  the  characteristic 
hump  in  the  data,  These  reduced  intensities,  scattered  and  unscattered, 
are  then  plotted  as  functions  of  beam  velocity.  The  plot  is  made  by 
computer  with  a  Calcomp  Model  593  plotter  which  produces  a  50  by  90  cm 
plot  with  precision  to  0.025  cm.  The  plot  is  used  in  a  smoothing 
technique  for  the  intensity  data  by  hand  drawing  the  best  smooth  curves 
through  the  scattered  and  unscattered  data  points.   The  coordinates  of 
several  arbitrary  positions  on  these  curves  are  then  obtained  using  an 
optical  centering  scope  on  a  verticle  milling  machine  which  has  table 
positioning  precision  to  0.25  nun.   This  technique  though  more  time  con- 
suming Droved  to  be  superior  to  the  several  computer  averaging  techniques 
which  vrere  attempted;  because,  unlike  these,  it  introduced  no  systematic 

bias  into  the  results. 
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These  smoothed  data  points  are  used  with  the  corrected  scattering 
gas  number  density  to  obtain  the  apparent  total  cross  section  as  a 
function  of  beam  velocity.  The  apparent  cross  sections  are  then  cor- 
rected for  motion  of  the  scattering  gas  by  using  the  Berkling  technique. 
The  resulting  values  of  this  correction  are  then  further  improved  by 
correcting  for  the  change  in  angular  resolution  due  to  the  change  in 
total  cross  section. 

A  plot  of  the  logarithm  of  these  cross  sections  as  a  function  of 
the  logarithm  of  beam  velocity  is  made,  and  the  best  portion  of  this 
is  used  in  determining  the  reactive  cross  section. 

Plot  Production 

Average  Transmitted  Beam  Velocity 

The  average  beam  velocity  during  an  intensity  measurement  is 

calculated  from  the  geometry  of  the  velocity  selector,  its  angular 

velocity,  and  the  temperature  and  mass  of  the  effusing  beam  using 

~  11 
the  following  integral 

where  FQJ)  is  the  transmission  function  which  is  dependent  on  the 
velocity  selector  geometry,  and  V.  and  V  are  calculated  limiting 
velocities  for  a  given  angular  velocity  of  the  selector. 

Intensity  Calculation 

Scattered  intensity  I  and  unscattered  intensity  I  are  calcu- 
lated in  the  same  fashion.  The  following  will  demonstrate  this  process 
using  just  IQ. 


25 
iQj)   =  P/[S  X(Wl 


where  F  is  a  normalization  factor  used  to  keep  the  values  of  IQ  between 
0  and  1,  S  is  a  scale  factor  to  account  for  different  output  ranges  of 
the  electrometer,  X  Qj)   is  the  charging  time  of  the  electrometer  at  the 
average  transmitted  velocity  V . 

In  order  to  increase  the  precision  of  the  plot,  the  intensity  is 
reduced  in  the  following  manner 


I0*(V)   =  I0QJ)  I     B5/2  exp  (- 


where 

,2 


=  my'  /   2  k  T, 


B 


for  which  m.  is  the  mass  of  the  beam  atoms,  k  is  the  Boltzmann  constant, 
B 

and  T-  is  the  temperature  of  the  beam  oven.  Figure  6  is  a  typical  plot 

B 

of  intensity  over  a  range  of  selector  velocities,  and  Figure  7  is  a 
typical  reduced  intensity  plot.   Values  of  I   and  Is  are  taken  at 
several  arbitrary  velocities  from  the  smooth  curves  drawn  through  the 
data.   So,  essentially,  the  plotting  technique  interpolates  the  data  to 
give  values  of  I  *  and  I  ^  at  the  same  velocity,  which  are  then  converted 
back  to  real  intensities  IQ  and  Ig. 

Scattering  Gas  Number  Density 

To  calculate  the  scattering  gas  number  density  distribution 
along  the  beam  axis,  the  number  density  of  the  source  of  gas  must  be 
known  at  a  point  near  the  slit-plugs  but  optically  occluded  from  the 
slits.   This  measured  number  density  is  then  corrected  for  nonuniformity 
along  the  beam  path  due  to  effusion  of  scattering  gas  through  the  beam 
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entrance  and  exit  slits.   The  results  of  these  calculations  appear  in 

the  expression  for  the  cross  section  as  a  correction  f actor  d .   It  has 

15,16,17 
been  shown        that  failure  to  introduce  this  correction  can  intro- 
duce an  error  of  several  percent  in  the  total  cross  section  and  even 
greater  error  in  the  reactive  cross  section.   Therefore,  its  inclusion 
in  this  study  is  essential,  and  since  its  calculation  is  new,  it  must 
be  described  in  detail.   That  description  is  given  in  the  appendix  of 
this  dissertation. 


Calculation  of  Source  Number  Density 

The  scattering  gas  pressure  is  monitored  with  a  thermocouple 
gauge,  which  shows  a  change  in  output  emf  with  a  change  in  pressure. 
The  method  for  determining  the  scattering  gas  pressure  is  quite  similar 
to  that  used  for  determining  intensity.   The  emf  data  are  plotted  versus 
beam  velocity  and  the  best  smooth  curves  are  drawn  through  the  data 
taken  when  the  box  was  evacuated  and  through  the  data  taken  when  the 
box  was  filled  with  scattering  gas.  The  ratio  R3Q/)  is  then  determined 
for  the  same  velocity  points  that  were  used  in  the  intensity  calculation 
by  the  equation 

reM  =  [EqCW  -  s(W]   /  20OV)    , 

where  EQQj)   is  the  emf  for  evacuated  conditions  and  EQJ)   is  the  emf 
for  filled  conditions.   Unfortunately,  the  relation  of  RE  to  pressure 
depends  on  the  nature  of  the  scattering  gas.   This  requires  that  for 
each  gas  the  thermocouple  gauge  be  calibrated  against  a  standard.   The 
calibration  is  carried  out  in  a  static  system  which  is  made  possible  by 
replacing  the  slit-plugs  in  the  scattering  box  with  sealed  plugs.  An 
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MICS  Baratron  Type  77  with  a  Type  77H-1  head  vras  used  as  a  standard. 
Pressure  P(l/)   is  found  by 

PQJ)  =  A  +  B  RE(W 

where  A  and  3  are  the  intecept  and  slope  respectively  of  the  straight 
line  obtained  by  least-squaring  the  pressure  calibration  data.   The 
number  density  then  is 

■  *0(U)   =  W)  I   k  Tg 

where  k  is  the  Boltzmann  constant  and  T  is  the  temperature  of  the 
scattering  gas. 

Total  Cross  Section 
Apparent  Total  Gross  Section 

The  apparent  total  cross  section  QA   is  obtained  from  the 
following  relation 

QApp(W  =  m  [i0(W  /  W)]  /a  nc  (l  +  t) 

where  L  is  the  distance  between  the  slit-plugs,  and  T  is  thickness  of 
each  slit. 

Berkling  Correction 

Because  the  target  molecules  are  not  at  rest,  the  apparent  cross 
section  will  be  too  large.   It  is  therefore  necessary  to  average  over 
the  distribution  of  relative  velocities  between  the  beam  and  target 
gas  molecules.   For  a  velocity  selected  beam  scattered  by  an  inverse 
power  attraction  by  a  Maxwellian  target  gas,  the  necessary  averaging 
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8 
has  been  computed  by  Berkling  et  al.   and  appears  in  the  expression 

for  the  cross  section  as  a  correction  factor  P  (s,  u/Uq)' 

lnllJV)  /   W)} 

Q  fM  = — — -7— 

e^f     cCnQ   (L  +  T)  F(s,  V/VQ) 

where  V     is  the  most  probable  velocity  of  the  scattering  gas.   Berkling 

G 
has  presented  a  complete  derivation  of  this  correction  factor  and  a 

rather  extensive  tabulation  of  F  values  for  various  values  of  s  and 

V/Vr.      Use  of  these  tables  greatly  simplifies  the  present  calculation. 

It  should  be  pointed  out  that  attempts  were  made  at  using  another 

18 
correction  scheme  worked  out  by  Kinsey,   but  the  results  were  unsatis- 
factory.  The  failure  arises  from  the  fact  that  the  Kinsey  method 
requires  a  very  precise  calculation  of  the  collision  relative  velocity, 
averaged  over  the  distributions  of  both  the  target  gas  and  the  selected 
beam  velocities.   Its  exact  calculation  involves  evaluation  of  the 
same  type  of  integrals  which  appear  in  the  Berkling  procedure.   Solution 
of  these  integrals  requires  extensive  use  of  a  high  speed  computer. 

If  one  neglects  the  narrow  soread  of  beam  velocities,  the  average 

19 
relative  velocity   is 

V  QJtVQ)   =  -I  (2  k  TQ  /  mG)  [2  exp(-ly  /  VQ   )  Iff' 
+  2  QJ/IQ   +  QJA/q)"1   erf  (Wg)]   , 

where  T  is  the  temperature  of  the  gas,  m  the  mass  of  the  gas,  and 
G  u 

erf  is  the  standard  Gaussian  error  function.   In  terms  of  this  relation 
the  cross  sections  were  calculated  for  Gs  -  Ar  by  both  the  Berkling  and 
the  Kinsey  methods.   The  results  (Figure  8)  show  a  divergence  between 
the  two  curves  at  low  velocities.  The  break  in  the  curve  for  the  Kinsey 
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corrected  cross  sections  resu] bs  from  the  finite  spread  of  beam  velocities 
Since  corrections  for  this  vrould  require  as  much  labor  as  computing  the 
Berkling  factors,  it  is  more  reasona.ble  to  use  the  3erkling  scheme  for 
the  present  work. 

Angular  Resolution 

This  last  correction  is  made  because  of  finite  angular  resolution 
of  the  apparatus.   This  arises  from  the  fact  that  the  scattering  region 
and  the  detector  wire  are  of  finite  size,  which  makes  it  possible  for 
particles  that  have  undergone  feeble  collisions  to  strike  the  detector 
and  be  counted  as  unscattered.   Tnis  is  especially  significant  when  the 

beam  ^articles  are  very  much  more  massive  than  the  target  particles  or 

7 
when  the  cross  section  is  large.   Kush  developed  a  scheme  for  deter- 
mining the  geometrical  resolution  0  of  the  machine  in  which 

G  =  wb  /  L  , 

where  W  is  the  half -height-width  of  the  beam  at  the  plane  of  the 
b 

detector  (hot-wire  filament)  and  L  is  the  distance  between  the  scattering 

20 
region  and  the  detector.   Kerschbach  et  al.   have  presented  all  the 

necessary  parameters  for  calculating  a  correction  factor  R(Q)  for  the 

apparent  cross  section  for  scattering  by  an  r   potential.  The 

resolution  correction  factor  is 

R(Q)  =  0.6248  -  0.4246  p  +  0.3752  exp(-  0.3129  p2) 

+  O.38I7  p  erf (0.9016  p), 

if  p  -S  1.592, 
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-1/? 
R(q)   -  0.0814  +  0.5604  p     /J 

if  p  ~  1 . 592  . 

The  resolution  parameter  p  is 

P  =  9m  ^Q*ff(W  /7rt  . 

where  m  is  the  mass  of  the  beam  particles.   As  defined,  the  resolution 
correction  factor  appears  in  the  denominator  of  the  cross  section 
equation 

In  I0(lO  /  W) 
^Ot0j)   =  F(s.U/V6)   R(Q)0CnQ  (L+  T) 

Since  the  correction  factor  is  dependent  on  the  cross  section,  and 
since  it  changes  the  value  of  the  cross  section,  the  process  of  cal- 
culating the  correction  factor  and  applying  it  to  the  cross  section 
is  iterated  until  convergence  of  less  than  0.01%  is  attained. 

Reactive  Gross  Section 

For  s  =  6,  taking  the  logarithm  of  equation  5  gives 

In  0^(1/)  =  -0.4  m^  +  3, 

where  B  is  a  constant.   Thus,  a  plot  of  In  %n(V)   verses  ln^ should 
give  a  straight  line  with  slope  equal  to 

-  0.4  =  -  2  /  (S  -  1). 

Plots  of  Qw^Cl/O  versus^*  for  reactive  systems  at  low  velocities 
have  slopes  less  than  -0.4.  The  reactive  cross  section  Q^  is  obtained 
by  reducing  Q,.,  QJ)   by  small  increments  and  finding  the  slopes  of  the 
resulting  logarithmic  plots  by  least- squaring  until  the  slope  is 
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exactly  -0.4.  At  this  point  the  sum  of  the  increments  of  Qm0t  is 
taken  to  be  the  reactive  cross  section  QD  . 

nX 


C^  '  Constant 


Using  the  intercept  of  the  final  logarithmic  plot  obtained  above, 
the  C    constant  is  calculated  by  the  following  equation 


C^6)  =  5.6?xlO"3°  exp(5/2  B)  , 


where  3  is  the  intercept  of  the  straight  line  obtained  by  the  least- 
squaring  Drocess.   The  units  for  the  constant  obtained  in  this  manner 

6 
are  erg  cm  , 


CHAPTER  IV 
RESULTS 

Discussion 

The  principal  purpose  of  this  research  project  was  to  deduce  a 

means  for  extracting  information  about  both  the  elastic  and  reactive 

channels  of  binary  molecular  interactions  fron  measurements  of  the 

velocity  dependence  of  the  absolute  total  scatteriag  cross  section, 

and  to  implement,  if  possible,  experimental  techniques  to  make  use 

of  the  separation  scheme.   A  completely  general  method  applicable  to 

all  possible  types  of  interactions  would  be  very  difficult  if  not 

impossible  to  construct.   Consideration  of  more  restrictive  special 

cases,  however,  holds  out  some  hope  for  at  least  limited  success.  The 

soecial  features  which  make  a  particular  interaction  suitable  for  this 

study  are   (l)  the  general  velocity  dependence  of  the  elastic  cross 

section  should  be  reasonably  predictable;   (2)  the  reactive  cross  section 

should  be  relatively  large;   (3)  the  probability  of  nonreactive  inelastic 

processes  should  be  low;  and  (4)  the  system  should  include  some  species 

that  is  easily  and  reliably  detected.  Fortunately,  it  happens  that 

these  features  are  more  or  less  simultaneously  met  for  a  large  class 

of  interactions  between  highly  reactive  gaseous  species.   As  early  as 

1 
the  1930* s  the  atomic  flame  experiments  of  Polanyi  showed  that  reactions 

between  alkali  metal  atoms  and  halogen  molecules  such  as 


K  +  Br2  =  KBr  +  Br 
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proceed  with  little  or  no  activation  energy  and  at  such  rapid  rates 

that  their  reactive  cross  sections  were  anticipated  to  be  in  excess  of 

100  1  ,  compared  to  the  10  %     or  less  of  nost  other  known  reactions. 

More  recent  ah  initio  calculations  of  potential  enerry  surfaces  for 

21 
such  ternary  systems  show  them  to  be  highly  attractive.    'The  energy 

barrier,  if  present  at  all,  is  low  and  occurs  early  in  the  course  of 

collision.   This  implies  that  the  reaction  probability  is  relatively 

high,  and  although  the  products  are  expected  to  be  formed  in  excited 

vibrational  states,  the  chance  of  vibrational  excitation  without 

reaction  is  small,  which  is  in  contrast  to  the  situation  resulting  when 

the  energy  barrier  occurs  late  in  the  course  of  reaction.  Finally,  the 

detection  of  alkali  metal  beams  by  surface  ionization  on  a  detector  of 

high  work  function  is  a  well  established  art,  so  this  general  type  of 

interaction  system  seems  uniquely  well  suited  for  the  purposes  of  the 

present  research.   Because  cesium  is  both  the  most  reactive  and  most 

easily  detected  alkali  it  was  chosen  for  the  beam  material.  Because, 

however,  the  elemental  halogens  are  quite  corrosive  and  therefore  hard 

on  the  apparatus,  some  less  objectionable  target  gases  were  desired. 

Other  research  has  shown  that  the  reactions  of  alkali  metal  beams  with 

halogen-containing  compounds  proceed  by  similar  mechanisms  to  those  for 

the  free  halogens.   Apparently,  aside  from  limitations  on  reactivity 

due  to  steric  hindrance,  the  particular  radical  to  which  the  halogen 

atom  is  bound  does  not  greatly  influence  its  extraction  by   the  alkali 

atom.   For  this  reason  it  was  decided  to  test  the  separation  scheme  on 

data  for  a  variety  of  halogenated  hydrocarbons.  The  plan  was  to  look 

for  trends  in  the  elastic  and  reactive  cross  sections  with  variation 

in  structural  features  (eg.  number  of  halogen  atoms,  kind  of  halogen 


atom,  isomeric  structure,  position  of  double  bonds,  etc.;.  Scattering 
measurements  for  several  systems  were  completed,  and  preliminary  analysis 
was  very  encouraging.  An  apparent  reactive  cross  section  was  detected 
for  each  system  in  which  reaction  was  expected.  A  nonreactive  system 
Cs  -  Ar  was  run  as  a  control  to  test  the  elastic  scattering  theory  and 
assure  that  the  apparent  reactive  cross  sections  were  not  artifacts  of 
the  experimental  procedure. 

Final  analysis  of  the  data,  including  estimation  of  the  limits  of 
uncertainty  in  the  results,  pointed  to  a  basic  flaw  in  the  experimental 
method  as  it  presently  exists.   The  extracted  value  of  the  reactive 
cross  section  is  very  sensitive  to  the  magnitude  of  the  total  cross 
section;  hence,  the  need  for  unusual  accuracy  in  the  measurements.  All 
of  the  innut  factors  in  the  equation  for  the  total  cross  section  are 
known  well  enough  to  introduce  a  net  uncertainty  of  less  than  \%  with 
the  exception  of  the  value  of  nQ,  the  number  density  of  scattering  gas 
in  the  source  ragion  of  the  scattering  chamber.   An  error  of  as  little 
as  2%   in  n~  can  lead  to  as  much  as  100?2  error  in  the  reactive  cross 
section.  For  most  of  the  experimental  runs  completed,  the  pressure 
data  from  which  n0  is  computed  showed  an  unusual  erratic  behavior.   This 
was  ultimately  traced  to  a  faulty  circuit  in  the  regulated  power  supply 
serving  the  thermocouple  gauge.   Even  with  the  circuit  repaired  the 
uncertainty  in  nQ   is  about  t  f%   This  could  be  further  reduced  to 
bring  this  measurement  more  in  line  with  the  other  parameters  by  appro- 
priate modifications  of  the  apparatus.   These  are  outlined  at  the  end 
of  this  section,  but  for  the  present  the  target  gas  pressure  measurement 
is  definitely  the  weak  link  in  the  cross  section  determination.   Because 
of  this  it  was  decided  to  present  here  only  the  most  reliable  data 
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presently  at  hand  to  indicate  the  potential  utility  of  the  proposed 
separation  method.   Recalling  that  order  of  magnitude  precision  is  all 
that  is  generally  available  for  reactive  cross  section  measurements, 
these  results  are  encouraging. 

The  system  which  was  most  studied  and  for  which  the  most  reliable 
data  exist  is  Cs-CILCICI^CI.   Three  separate  runs  on  this  system  were 
made.  The  possible  precision  of  the  method  suffered  from  imprecision 
in  the  pressure  measurements. 

C/  '   Constant 

Table  2  gives  values  of  the  CT  ;  constant  obtained  from  three 
separate  measurements  of  the  Cs-CH^CICH-Cl  interaction  using  three 

techniques  of  smoothing  the  scattering  gas  pressure  data.  The  overall 

-57       6 
average  value  obtained  from  these  measurements  is  2. 34x10  ■"    erg  cm 

with  an  average  deviation  of  ±.9%.     For  the  nonreactive  system  Gs-Ar 

the  C  ^constant  obtained  is  5.4SxlO~^  erg  cm6  with  an  average 

deviation  oftz%.      It  should  be  noted  that  the  Qr° '   constant  obtained 

for  this  system  without  the  use  of  the  ot  correction  factor  is  4.60x10 

erg  cm.  Thus,  in  this  case  explicit  accounting  for  the  nonideality  of 

the  target  gas  number  density  changes  the  value  of  the  Van  der  tfaals 

constant  by  16%,     an  amount  which  if  not  corrected  could  change  the 

calculated  value  of  Q   by  several  orders  of  magnitude. 
Rx 

"  Total  Gross  Section 

Table  3  gives  values  of  each  factor  used  in  equation  7  to 
calculate  the  absolute  total  cross  section  at  several  velocities  for 
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TABLE  2 


Experimental  Values  of  the  G^  '   Constant  for  Cs-CHgClGh^Gl . 


Scattering  Gas  Pressure 
Data  Smoothing  Technique 


Experimental  Hun  Number 


I663      1691      1693 


Hand-fit  curve 

Least- squared  straight  line 

Averaged  constant  value 


2.13 

2.59 

2.56 

2.20 

2.4o 

2.51 

2.03 

6  (  <rr57\ 

All  values  are  in  erg  cm   (.xlO   ; 
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TABLE  3 


Typical  Velocity  Dependent   Intermediate  Results  and  Parameters 
Leading  to  Total   Cross  Sections, 


Selector  Velocity  RPM  2000  3000  4000  5000  6000  7000 

dn0(L  +  T)  cm-2(xl012)  1.930  1.920  1.925  1.928  1.930  1.935 

Beam  Velocity       cm  sec^UlO^)     1.457  2.179  2.893  3. 597  4.269  4.968 

QApp(W  &  3335  2399  1922  I636  144b  1302 

F(s,  V/VG)  1.4778  1.2530  1.1509  1.0981  I.0687  I.0509 

%ervS^  f'2  225?  191k  16?°  lZ+9°  1355  1?'39 

r(q)  0.9222  0.6926  0.8672  0.8442  0.8232  0.6046 

Qrot(19)  A2  2447  2144  1926  1765  1646  151+0 


44 

one  experiment.  Also  shown  are  intermediate  results  of  the  application 
of  individual  correction  factors.   It  is  interesting  to  note  that 
application  of  the  Berkling  correction  factor  has  reduced  the  cross 
sections  and  that  application  of  the  angular  resolution  correction 
has  increased  them. 

Reactive  Cross  Section 

Table  4  shows  values  of  Q   for  Cs-CHgClCHgCl  for  three  separate 
experiments  using  several  techniques  for  smoothing  the  scattering  gas 
pressure  data.  The  average  value  of  Q^  obtained  using  all  these  data 
is  115  1  6°  £2. 

From  this  value  of  the  reactive  cross  section  the  reaction  rate 

14   3   ,  -1    -1 
constant  calculated  from  equation  5  is  about  2  x  10   cm  mole   sec   . 

1 
Using  the  "harpoon  mechanism"  first  envisioned  by  Polanyi,   the 

electron  affinity  of  C^CICM-Cl  can  be  approximated  uy  the   relationship 


2 
e 


R  - 
x 


(I. P.  -  S.A. ) 


where  I. P.  is  the  ionization  potential  of  cesiura,  and  S.A.  is  the 
electron  affinity  of  CHg&CHgCl.   Rx  is  the  larSest  separation  where 
an  electron  can  transfer  from  the  alkali  metal  atom  to  the  halogen- 
containing  molecule.   Since  one  of  the  features  of  the  harpoon  scheme 
is  that  reaction  immediately  follows  this  charge  transfer,  it  is 
assumed  that 


2 

*Rx     x 


TABLE  it- 


Results  of  Q   for  Cs-GHoClCHpGl  Determinations  for  Three  Separate 
Experiments  losing  Various  Scattering  Gas  Pressure  Smoothing  Techniques. 


Scatteriag  Gas  Pressure 
Averaging  Technique 


Hand-fit  curve 

Least- squared  straight  line 

Averaged  constant  pressure 

Hand-fit  curve  over 
extended  velocity  range 

Least- squared  straight  line 
over  extended  velocity  range 


Experimental  Run  Number 


I663 


1691 


120 


1693 


90       24       105 

54       51 

298      48       310 


79 


9C 


All  values  of  Q.?^  are  in  K   . 


i+6 


The  value  of  the  electron  affinity  of  CK  CI  CI'  CI  calculated  using  the 

present  results  is  1,3  ev. 

For  the  purpose  of  demonstrating  where  the  weak  point  in  the 

experiment  lies,  consider  the  entries  in  table  5*   It  is  obvious  that 

better  measurements  of  scattering  gas  pressure  are  needed  to  improve 

Q   determinations.  3y  a  propagation  of  errors  analysis  it  was 
Rx 

discovered,  as  was  mentioned  earlier,  that  as  little  as  a  2%   change  in 
scattering  gas  number  density  will  cause  a  10Q%  change  in  the  calculated 
reactive  cross  section  for  this  system. 

Discussion  on  Improving  Pressure  Measurements 

The  pressure  measurement  procedure  involves  writing  down  values 

of  the  thermocouple  gause  emf  E  and  E  for  each  measurement  of  unscat- 

0 
tered  and  scattered  beam  intensity.   One  must  monitor  the  emf  manually 

with  the  potentiometer  circuit.   Small  fluctuations  in  the  output 

signal  during  steady  state  conditions  prevent  the  operator  from 

deciding  unambiguously  on  the  magnitude  of  the  output  signal. 

The  pressure  and  hence  n„  is  related  directly  to  the  fractional 

change  in  emf  on  introduction  of  scattering  gas  fvE  -  S)  /  S  =  Rg. 

Because  the  pressures  required  to  maintain  the  gas  density  at  levels 

appropiate  for  scattering  20  -  2%  of  the  incident  beam  are  low 

(ca.  10~^  torr),  the  corresponding  values  of  R  are  low  (ca.  0,5%)', 

that  is,  R  depends  on  a  factor  which  is  the  difference  between  two 

numbers  which  are  nearly  equal.   For  the  actual  gauge  used,  the  value 

of  E  is  about  18  mv.   An  uncertainty  in  the  values  of  Eq  and  S  of 

as  low  as  t   10  volt  can  thus  introduce  a  2%  uncertainty  in  the  value 

of  n  .   This  problem  is  somewhat  like  that  which  once  existed  in  the 
0 
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TABLE  5 
Precision  of  Input  Data 


Intensity  G.001; 

Beam  oven  temperature  '  0.1% 

Scattering  gas  pressure  l/o 

Scattering  box  temperature  0.1% 
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measurement  of  beam  intensities.   tfnen  IQ  and  I  are  observed  directly, 
there  are  always  snail  fluctuations  which  introduce  uncertainty  in 
direct  determination  of  their  magnitudes.   This  difficulty  vras  overcome 
by  changing  to  an  indirect  node  in  which  the  signal  is  averaged  over 
a  period  of  time  (ca.  5  sec)  that  is  short  with  respect  to  the  time 
scale  of  the  experiment  but  long  compared  to  the  fluctuation  frequency. 
Using  this  method  with  the  solid-state  circuitry  described  in  the 
exuermental  section,  it  is  possible  to  measure  beam  intensities  to  a 
precision  of  better  than  one  part  in  1(P.   The  implication  is  clear 
that  a  similar  integrating  procedure  applied  to  the  thermocouple  gauge 
output  could  be  expected  to  upgrade  the  precision  of  that  measurement 
significantly.   It  is  recommended  that  such  a  modification  to  the 
apparatus  be  made  so  as  to  permit  full  extraction  of  its  potential 
utility. 


APPENDIX 
CALCULATION  OF  & 

The  loss  in  intensity  of  a  monoenergetic  beam  scattered  in  a 
homogeneous  gas  of  number  density  n  and  path  length  L  is  given  by  the 
familar  Beer- Lambert  lav:  as 

ln(lQ  /  I)  =  Q  nQ  L  , 

where  Q  is  the  total  scattering  cross  section.  When  the  scattering  gas 
number  density  is  not  uniform  along  the  beam  path, this  equation  must  be 
written 


ln(l0/l)  =  QfDn 


zdZ  ' 

where  n„  is  the  target  gas  density  near  point (Z) centered  on  the  beam 
axis,  and  the  integral  is  evaluated  along  the  beam  (z)  axis  from  the 
position  of  the  source  to  that  of  the  detector.   When  the  scattering 
gas  is  introduced  to  a  scattering  chamber  located  somewhere  along  the 
beam  path,  its  steady  state  distribution  is  always  nonuniform  due  to 
effusion  through  the  beam  passage  slits.   In  the  ideal  case  for  which 
the  slit  walls  are  assumed  to  radiate  gas  at  a  level  characteristic  of 
walls  in  equilibrium  with  gas  with  number  density  nQ,  it  has  been  shown 
that  the  value  of  the  integral  approaches  the  value  L  +  T,  where  T  is 
the  slit  channel  thickness.    Because  in  actual  scattering  chambers  the 
effective  density  deviates  from  the  ideal  case,  it  is  necessary  to 
account  for  this  effect  for  accuracy.  This  is  done  by  computing  the 
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actual  integral  after  the  steady  state  density  distribution  has  been 
determined.  To  simplify  the  final  notation,  this  integral  is  expressed 
as 


nr  dZ  =cC(L  +  T)   , 


from  which  one  finds 


f 


L  +  T 


The  procedure  for  calculating  cC  for  the  scattering  chamber  used  in  this 
work  is  given  below. 

The  calculation  of  the  density  distribution  correction  factor  <£ 
for  the  scattering  chamber  used  in  this  work  (Figure  2;  uses  the  results 
for  the  semi-infinite  geometry   as  a  basis,  with  finite  slit  height  and 
finite  channel  thickness  treated  as  perturbations.  Figure  9  is  a  schema- 
tic representation  of  the  scattering  chamber  of  finite  slit  height  and 

thickness . 

16 

Following  the  procedure  outlined  by  Colgate,   the  faces  of  the 
slit-plugs  are  divided  into  uniform  parallel  bands,  i  on  one  and  j  on 
the  other.  All  lengths  are  reduced  in  terms  of  the  slit  width  (i.e. 
the  slit  width  is  taken  as  one  unit).   The  source  region  is  considered 
to  have  a  uniform  gas  number  density  nQ  =  1  throughout.   The  bands  are 
of  width 

Ax=(S-  0.5)  /  m  , 

where  m  is  the  number  of  bands  chosen.   Each  band  is  assumed  to  have 
uniform  effective  radiative  density  ry   Initial  values  of  ni  assigned 
to  each  band  are  the  results  for  the  semi-infinite  case.  Details  for  the 


Figure  9.   Schematic  of  scattering  box  slit  geometry. 
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calculation  of  the  effects  of  finite  channel  thickness  and  slit  height 
follow. 

Finite  Slit  Height 

The  major  effect  of  this  change  is  that  caused  by  the  new  areas 
of  source  above  and  below  the  slit.  The  contribution  to  the  flux  of 
molecules  onto  the  center  of  the  jth  band  resulting  from  the  additional 
radiation  beyond  y  =  t  H  from  the  ith  band  on  the  opposite  wall  is 


A 


-..=  (1-n.)  Q  (eii-isin2eij) 

hrVC 


xi.i  +  ^x 


(xlj  +|Ax)2  +  L2 


x.  .  -  -§Zx 


(xi<   -  i^x)     +  L 


where 


-V 


ij 


wztttt? 


ij 


and 


x.  .  =     x.    -  x. 
ij  i         0 


4x 


for  bands  on  the  same  side  of  the  slits,  and 


x  .  =  x.  +  x.  -1 1  Ax  -  1 

ij   J  i    J   | 

for  bands  on  opposite  sides  of  the  slits.   The  term  V   is  the  average 
velocity  of  the  gas  particles.  The  contribution  to  the  flux  onto  the 
center  of  the  jth  band  due  to  source  above  and  below  the  slit  is  given  by 

A  IT 4  =£_  (ft  -  t  sin2ei)  (Asm©.)  , 


r/> 


where 


and 


~  -1 

0,  =  cos   Ur?-      ?       z" 

J  l  VH    +  L    +  x 

j 


x.  +  t  x. 

A  sin0.  =     ,        -1  -7=^ 


j   Vu7TT)T7  V(x,  -  i)  +  ^ 


J  J 

Using  the  results  of  these  calculations,  new  effective  densities  are 
calculated  for  each  band.  The  change  in  the  number  density  is  given  by 

l   "*■    toJ     -lj 

0 

and  the  new  radiative  density  is 

o   . 
n.  =  n  +&  n  , 
i    i     i 

where  the  n.  's  are  the  initial  trial  radiative  densities.   These  new 

l 

n.'s  are  then  used  as  trial  values  to  recalculate  theAn.'s,  and  the 
l  i 

process  is  iterated  until  satisfactory  convergence  is  attained. 

Finite  Channel  Thickness 

The  approach  to  this  perturbation  is  similar  to  that  used  above. 

Calculate  the  effective  radiative  density  on  the  channel  walls  due  to 

contributions  from  the  source,  scattering  chamber  walls,  and  the  channel 

walls  themselves.  Then  determine  the  effects  of  radiation  from  the 

channel  walls  onto  the  scattering  chamber  walls,  and  continue  by  iteration 

to  convergence. 

The  contribution  n    due  to  source  is 
CHS 
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L  +  T/2 


T/2 


JHS 


L+  T/2)2  +   (S  +Y)2f       (i  T2  +  1)2 


+  1/77- (6.  -  I  sin20_)     1 


L  +  T/2 


|(L  +  T/2)2  +  1 


where 


0     =  cos 


-1 


*[h2  +  (L+  T/2)2  +  i 


The  contribution  n^,tTTT  due  to  radiation  from  the  scattering  chamber 
CHW 


walls  is 


where 


nnuu  =  1/rrT,  n.    (^2  -0.  +1  sin20. )  A(sin0)     , 
CHW  /tii  i  i  i 


-1 


ei  =  C0S     \  [h2  +   (L  +  t/2)2  +   (x.    +  i)2] 


and 


A  (sin6).   = 


L  +  T/2 


L  +  T/2 


[(L  +  T/2)     +   (x.    -  |Ax)2]  2  [(L  +  T/2)2  +   (x.   +  |A 


x)' 


21* 


The  contribution  n„ r„  due  to  radiation  from  the  other  channel  walls  is 
CriC 

n         =  1/fr    (n         +  n        +  n       )  {  (*/2  -  9PH1   +  |  sin2£       ) 


GHG 


GKS  CHW  CHC 

T 


Xtt/h  +  l)2 


+   (^2  -6PHO  +1  sin2a„9) 


L  +  3/2  T 


L  +  T/2 


[(L  +  3/2  T)2  +  l]  *  [(L  +  T/2)     +  l]  2/j 


'CHI 


;32  +  D1 


and 


'CH2 


3  +  (L  +  T)  +  1  2 


Thus,  the  effective  radiative  density  n   is  the  sun  of  the  above 
contributing  terms,  or 


nGH  =  "cHS  +  "CHW  +  nCHC 


Now  the  effect  of  radiation  from  the  channel  walls  onto  the 
scattering  chamber  walls  is  determined.   The  change  in  density  in  the 
ith  band  is  given  by 


where 


and  for  which 


AnGH,i  =  nGH  Gi  ' 


1   n  oi         1 


Fei  =TT/2  -0i  +  1  sinaOi  , 


-1 


['H2+  (x.+|)2+  (L  +  T/2)2]*/ 


and 
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A  sir   p    ~ 


X.     +  ■';■ 
1 


X      + 
1 


i  "     [(xi+t)2  +  L2]1         [(x1+|)2+   (L  +  T)2]1 
To  simplify  the  iteration  process  the  ratio 


R.    = ^£ 

n 

is  calculated  for  each  band.   New  values  are  then  calculated  for  this 
ratio  by  the  following: 


6  n.  .   R.[g.  .  (Asin0).  .1 


where 


(Asintf 


CH 


X.  .  +  yAx 
_i3_ 


2  ,  TZ 


ij     [(x^  +  iaxr  +  i; 


X    -J^X 

— iJ . — 


(x.  .  -  iAx)  +  L' 


and 


for 

which 

0.  .  =  cos 

IJ 

The 

new 

value 

of 

R 
i 

is  thus, 

G.  .  =Tr-0.  .  +  |  sin20.    , 


H 


L  (H2  +  L2  +  x  2  )2  I 
ij 


R  =  C  +  2i^n.  .  /  n  ) 
i    l   j"  ij'  Chi 


These  new  values  of  R^  are  then  used  for  calculating  new  values  of 


An   /  n  ,  and  the  iteration  process  is  continued  until  satisfactory 
ij  '   CH 


58 

convergence  is  achieved.   With  the  final  values  of  R  new  values  of  the 
cell  densities  on  the  scattering  chamber  walls  are  calculated  by  the 
equation 

n.  =  n.°  +  n   R.   , 

ll      Ln   i 

where  the  n,°'s  are  the  previous  values  of  the  ru's.   These  n/s  are 
then  used  to  calculate  a  new  value  of  nnu,   which  is  then  used  in  the 
equation  above  to  recalculate  the  n^s.  This  process  is  repeated  until 
the  desired  convergence  is  reached.   Each  iteration  is  terminated  when 
agreement  between  successive  values  is  better  than  one  part  in  10  . 

Number  Density  Along  Beam  Axis 

16 
Using  the  equation  set  forth  by  Colgate, 

np  =  _1_    / nq  dw, 
k-'f  J    b 

where  n  is  local  number  density  near  point  P  on  the  beam  axis,  n  is 
P  a 

the  effective  radiative  density  at  surface  element  dS,  and  dw  is  the 
solid  angle  subtended  by  dS  at  P;  the  number  density  distribution  along 
the  entire  beam  path  can  be  calculated.  For  this  case  nQ   is  the  above 
calculated  value  n. ,  and  dw  the  associated  solid  angle.   For  simplicity 
of  calculation  the  process  is  divided  into  two  steps.   The  first  is  for 
contributions  to  number  density  inside  the  scattering  box,  and  the 
second  is  for  contributions  to  number  density  outside  the  box. 

Inside  Number  Density 

The  z-axis  inside  the  scattering  chamber  is  divided  into  K  equal 
segments  of  length  A?  =   M  /  L.   With  M  chosen  as  an  even  number,  the  jth 
segment  is  centered  at 
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Pj  =  (j  -  i)AP. 

T^e  contribution  to  density  at  P..  due  to  radiation  from  the  source  in 
the  x-diroctiona  ia  £iven  by 


n  (J)  =  lAf 

SI 


-1       1 
tan"1 (P-)  +  tan  X(L  -  P.)   . 


-a  n"**  I 


The  contribution  from  source  above  and  below  the  slit  (i.e.  in  the  y- 
directions  ard  with  lyl  =*  H)  but  not  in  line  with  the  slit  is 


,->   <f(wS2(i'j>' 


S2 


where 


wqo(i,j)  =^ll 


'S2 


o  o  2 

I   (JT  +  P/  +  x.  ) 


,1 


,  xi  +|Ax 

tan  —  -  tan 


-1  Xi  -^x 


The  contribution  to  number  density  at  P.  due  to  source  in  line  with  the 
slit  and  for  lyl  *■  H,  again  for  the  P  =  0  plug-face,  is 


n  (j)  =11- 
°>  ft 


L(H  +  p.2)2J 
J 


tan"1(lL 

■  2  P, 


The  contribution  at  P  due  to  radiation  from  the  ith  band  on  the  plug- 

J 
face  at  P  =  0  in 


B     £tT  i=l  x  b 


where 


KB(i,J)  =  <* 


,,2     2     ?J 
I (H  +  P.  +  x   ' 
0     i 


tan  1  i 


1  **  +  ^11  -  tan-^i'^") 
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Finally,  the  contribution  to  number  density  at  P,  due  to  the  wall  surface 
which  is  in  line  with  the  slit,  and  above  and  below  it,  for  the  plug-face 
at  P  =  0  is  given  by 


"1 

4     ^ 


.  (r  +  p.  Y     (b   +  i 


3    ''2 


-Z-X  tan 


2  P. 
J 


where  n  is  the  number  density  in  the  first  band  next  to  the  slit.  Thus, 

1 
the  total  axial  number  density  for  the  region  inside  the  scattering  box  is 


*5 


•1/2 


Inside     j=l 


A  L  nsi 


nq1(j)  +  n__(j)  +  n  (l-H-l-j)  +  n  (j)  +  n  (H+1-; 


S2 


S2 


S3 


3) 


B 

Outside  Number  Density 


.+  n  (j)  +  n  O'H-l-j)  +  n  (j)  +  n,(l'H-l-j)]AP  . 


The  calculation  of  the  number  density  distribution  outside  the 
scattering  box  is  carried  out  in  a  similar  fashion.   For  simplicity  the 
calculation  is  divided  into  two  parts;  one  for  P  's  that  have  direct 
radiative  contribution  from  the  source,  and  one  for  P.'s  that  have  no 
direct  contribution  from  the  source.  The  limits  are  defined  by 


I  (?.   -  L) 


where 


x.  +  |  A  x 

P.  = 

3       x.  +|Ax  -  | 


If  Q.  =»  H,  some  source  is  visible  to  P.,  and  if  Q.  =  H,  no  source  is 
3  3  3 

visible  to  P  .  For  this  outside  calculation  there  is  one  P  for  each 
3  3 
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ith  band  on  the  scattering  box  wall,  and  P.  is  the  most  distant  from  the 


slit. 


For  Q  =-  H  the  contribution  to  number  density  at  P.  due  to  radia- 


tion from  the  source  is 


where 


n  (j)  =  w  (j)   , 
S     __$ 


w  (j)  sr  4  tan 


B  /  2 


(P.  -  L)  |  +  B  +  (P.  -  if 


-  h   tan" 


(x.  +  | Ax)  H 


Pj  [(x.  +|/\x)2  +  H2  +  P^2]' 


The  contribution  due  to  radiation  from  the  wall  bands  is 


-.")  ■  £ 


n,  wT)(i,j) 


i B 


^rr 


4Tf 


^  tan 


-1 


H  /  2 


Pj  (i  +  h2  +  Pj  y 


-  4  tan 


-1/ 


3/2 


'i  +   B2  +  P.2)'/ 


where 


w„(i,j)  =  ^  tan 


-1 


(x.  +  |Ax)  H 


P.  [(xi  +  |Ax 


i  A  Y>2  ,  u2  ,  P  2)'i 


-  k   tan 


(x.  -Ux)  H 


P  |(x.  -  i- Ax)2  +  H2  +  P.2 

1  j  L  l  J 
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For  0  :f  H,  the  contribution  from  the  source  n  (j)  =  0,  and  the  contribution 

S 
from  the  bands  is  given  by 


nB(j; 


a  rn-  w  Cit j)i   ni  f       < 


Q.  /  2 

. J_ 

p.  (t  +  q/-  +  p.  y 


4  tan 


-1 


B  /  2      \1 


2     2  i 

p.  (t  +  3  +  Fj  r 


Thus , 


nTot(j)  =  nS0)  +  nBG)  ' 


and  the  total  number  density  outside  is 


n       =2 

Outside 


P  -  L 
m 


m 

n_  .,  ■  (1)  +  £  I(j)   , 
inside       M 


v;here 


nT  .„  (D=2  fnq1(l)  +  n  (1)  +  n  (M)  +  n  (1 )  +  n  (M)  +  n  (1)  +  njil) 


Inside 


SI 


S2V  '  SZy    '  S3      s3 


+  \^) +  VM)J  • 


The  l(j)*s  are  given  by 


and 


I(o)=nTot(o)(^L^l  .   «J-1. 


IP  -1  —  P •    \ 
-J £i  ,  if  j  ^  i. 
2     ' 


The  number  density  beyond  this  is  given  by 


63 


2nivL  -i  3/2 


tan        ,_  _^fZ  .^    .    „2 


Tail"   Tf    YL    "     \(P    -  L)V(P.  -  l)    +r/  +  t 

i  i 


'iVT2  +  B^  +  i  /J 


where  A P  =  1  and  P.    =  P.    +  iAP  . 
l  i 

Calculation  of  d    ■ 

Since  it  is  desired  to  express  the  final  result  in  terns  of  the 
correction  factor  oL,    the  separate  results  above  are  combined  as  follows 

*  =  (-inside  +  "outside  +  Vil  +  2   "cH  ^   (L  +  T> 

Table  6  presents  values  of  d   for  ideal,  knife-edge  channel,  and 
finite  channel  cases  for  the  particular  geometry  used  in  this  experiment, 
while  Figure  10  shows  the  number  density  profile  along  the  beam  axis  for 
the  scattering  box  used  in  this  experiment. 
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TABLE  6 


Values  of  n's  and  d.   for  Static,  Knife-edge  Slit,  and  Finite  Channel  Slit 
Cases  for  the  Scattering  Box. 


Inside 


Outside 


Taxi 

Channel 
ec- 


static     iCnife-edge   Finite-channel 
slit  •        slit 


8. 036 

6.695 

7.064 

0.9^9 

0.769 

0.810 

0.267 

0.231 

0.235 

0.500 

0.00 

0.458 

1.001 (1) 

a      (2) 
0.O55      ' 

('1  ) 
0.877v   ' 

(1)  L  +  T  =  9-66 

(2)  T  =  0;  thus,  L  +  T  =  L  =  9.26. 
All  n's  are  reduced  such  that  nQ  =  1. 


Figure  10,   Profile  of  scattering  gas  number  density  along  the  beam  path,. 


(,G 


P(J) 


BIBLIOGRAPHY 

1.  M.  Polanyi,  Atonic  Reactions,  Williams  and  Northgate  Publishing 
Company ,  London  (1932 J . 

2.  J.I.  Steinield  and  J.  L.  Kinsey,  in  Progress  in  Reaction  Kinetics, 
vol.  5,  G.  Porter,  alitor,  Pergamon  Press,  London  (I9?0), 

3.  J.  I.  Steinfeld  and  J.  L.  Kinsey,  ibid. 

4.  I.  Amdur  and  G.  Hammes,  Chemical  Kinetics;  Principles  and  Selected 
Topics,  McGraw-Hill,  New  York  (1966"]", 

5.  J.  B,  Hasted,  Physics  of  Atomic  Collisions,  American  Elsevier  Co., 
New  York  (1972). 

6.  R.  D.  Levine  and  R.  3.  Bernstein,  Molecular  Reaction  Dynamics, 
Oxford  university  Press,  Hew  York  (1974). 

7.  P.  Kush,  J.  Chem.  Phys . ,  40,  1  (1964). 

8.  K.  Berkling,  R.  Helbing,  K.  Kramer,  H.  Pauiy,  C.  Schlier,  and  P. 
Toschek,  Z.  Physik,  166,  406  (1962). 

9.  L.  D.  Landau  and  S.  H.  Lifshitz,  Quantum  Mechanics,  Pergamon  Press, 
London  (1959). 

10.  R.  C.  Schoonmaker,  J.  Phys.  Chem.,  6j,  892  (1961). 

11.  S.  0.  Colgate  and  T.  C.  Imeson,  Rev.  Sci.  Inst.,  3_6»  932  (1965). 

12.  S.  0.  Colgate  and  T.  C.  Imeson,  J.  Chem.  Phys.,  5J,  1270  (1970). 

13.  T.  C.  Imeson,  Ph.  D.  dissertation,  University  of  Florida  (1965). 

14.  C.  B.  Smith,  Ph.  D.  dissertation,  University  of  Florida  (1970). 

15.  S.  0.  Colgate,  Vacuum,  21,  483  (1971). 

16.  R.  N.  Nelson  and  S.  0.  Colgate,  Phys.  Rev.,  8,  3045  (1973;. 

17.  3.  P.  Hathur,  J.  S.  Field,  and  3.  0.  Colgate,  Phys .  Rev.,  11  (1975). 

18.  N.  C.  Lang,  H.  V.  Lilenfeld,  and  J.  L.  Kinsey,  J.  Chem.  Phys . ,  5J5., 
3114  (1971). 

19.  E.  H.  Kennard,  Kinetic  Thepry_  of  Gases ,  McGraw-Hill,  New  York  (1936). 


67 


20.  R.    J.    Cross,    S.   A.    Gislason,   and  D.    R.    Herschbach,    J.    Chem.    Phys . , 
45,   3582   (1966). 

21.  C.   F.    Bender,    P.    K.    Pearson,    S.    V.    O'Neil,    and  H.    F.    Schaefer,    J. 
Chen.    Phys.,    $6,   4626    (1972). 


BIOGRAPHICAL  SKETCH 


Jay  E.  Field  was  born  on  November  4,  19^7 »  in  Roanoke,  Virginia. 
He  was  brought  up  and  graduated  from  high  school  in  Norfolk,  Virginia. 
In  June,  1970,  he  received  a  Bachelor  of  Science  degree  in  Chemistry  at 
Virginia  Polytechnic  Institute  in  Blacksburg,  Virginia.   In  September, 
1970,  he  entered  graduate  school  at  the  University  of  Florida,  where  he 
worked  as  a  graduate  teaching  assistant  until  August,  1975* 

He  is  married  to  the  former  Laura  Lee  Alley,  and  they  have  a 
daughter,  Amy  lee. 


69 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


'/Z?W! 


f  /)  /Jr,L 


Samuel  0.  Colgate,  Chairman  / 
Associate  Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


n 


X, 


Marvin  "L.  iluga 

Associate  Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Fnilosophy. 


Xvn   £-•  ^& 


vjO.  3.  \±, 


Willis  3.  Person 
Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  ny  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Pnilcsophy. 


£Wib 


JeoVge  E.  Ryschkewitsch 
Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


John  W.  Flowers 

Professor  of  Physics  and  Astronomy 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  Department 
of  Chemistry  in  the  College  of  Arts  and  Sciences  and  to  the  Graduate 
Council,  and  was  accepted  as  partial  fulfillment  of  the  requirements  lor 
the  degree  of  Doctor  of  Philosophy. 

August,  1975 


Dean,  Graduate  School 


UNIVERSITY  OF   FLORIDA 


3  1262  08553  3205 


